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ABSTRACT 

We calculate the linear torque exerted by a planet on a circular orbit on a disc con- 
taining a toroidal magnetic field. All fluid perturbations are singular at the so-called 
magnetic resonances, where the Doppler shifted frequency of the perturbation matches 
that of a slow MHD wave propagating along the field line. These lie on both sides of 
the corotation radius. Waves propagate outside the Lindblad resonances, and also in 
a restricted region around the magnetic resonances. 
, The magnetic resonances contribute to a significant global torque which, like the 

Lindblad torque, is negative (positive) inside (outside) the planet's orbit. Since these 
resonances are closer to the planet than the Lindblad resonances, the torque they 
CO . contribute dominates over the Lindblad torque if the magnetic field is large enough. In 

■ addition, if (3 = c 2 jv\ increases fast enough with radius, the outer magnetic resonance 

becomes less important and the total torque is then negative, dominated by the inner 
O ■ magnetic resonance. This leads to outward migration of the planet. Even for (3 ~ 100 

I 1 at corotation, a negative torque may be obtained. A planet migrating inward through 

a nonmagnetized region of a disc would then stall when reaching a magnetized region. 
It would then be able to grow to become a terrestrial planet or the core of a giant 
planet. In a turbulent magnetized disc in which the large scale field structure changes 
sufficiently slowly, a planet may alternate between inward and outward migration, 
depending on the gradients of the field encountered. Its migration could then become 
i diffusive, or be limited only to small scales. 

5-j \ Key words: accretion, accretion discs - MHD - waves - planetary systems: proto- 

planetary discs 



1 INTRODUCTION 

Almost 20% of the extrasolar planets detected so far orbit at a distance between 0.038 and 0.1 astronomical unit (au) from their 
host star. It is very unlikely that these short-period giant planets, also called 'hot Jupiters', have formed in situ (Bodenheimer 
1998; Bodenheimer, Hubickyj & Lissauer 1999). More likely, they have formed further away in the protoplanetary nebula and 
have migrated down to small orbital distances. It is also possible that migration and formation were concurrent (Papaloizou 
& Terquem 1999). 

So far, three mechanisms have been proposed to explain the location of planets at very short orbital distances. One 
of them relies on the gravitational interaction between two or more giant planets, which may lead to orbit crossing and to 
the ejection of one planet while the other is left in a smaller orbit (Rasio & Ford 1996; Weidenschilling & Marzari 1996). 
However, this mechanism cannot account for the relatively large number of hot Jupiters observed. Another mechanism is the 
so-called 'migration instability' (Murray et al. 1998; Malhotra 1993). It involves resonant interactions between the planet and 
planetesimals located inside its orbit which lead to the ejection of a fraction of them while simultaneously causing the planet 
to migrate inward. Although such interactions may lead to the migration of giant planets on small scales, they would require 
a very massive disc to move a Jupiter mass planet from several astronomical units down to very small radii. Such a massive 
disc is unlikely and furthermore it would be only marginally gravitationally stable. The third mechanism, that we are going 
to focus on here, involves the tidal interaction between the protoplanet and the gas in the surrounding protoplanetary nebula 
(Goldreich & Tremaine 1979, 1980; Lin & Papaloizou 1979, 1993 and references therein; Papaloizou & Lin 1984; Ward 1986, 
1997). Here again the protoplanet can move significantly only if there is at least a comparable mass of gas within a radius 
comparable to that of its orbit. However this is not a problem since this amount of gas is needed anyway in the first place to 
form the planet. 
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Note that hot Jupiters can also be produced by the dynamical relaxation of a population of planets on inclined orbits, 
formed through gravitational instabilities of a circumstellar envelope or a thick disc (Papaloizou & Terquem 2001). However, 
if objects as heavy as r-Boo may be produced via fragmentation, it is unlikely that lower mass objects would form that way. 

A planet embedded in a gaseous disc launches density waves at the Lindblad resonances (Goldreich & Tremaine 1979, 
hereafter GT79). The torque exerted by the planet on these waves is responsible for the exchange of angular momentum 
between the disc rotation and the planet's orbital motion. The interaction between the planet and the disc inside/outside its 
orbit leads to a negative/positive torque on the disc, and therefore to a gain/loss of angular momentum for the planet. In the 
linear regime, because in a (even uniform) Keplerian disc the outer Lindblad resonances are slightly closer to the planet than 
the inner Lindblad resonances, the interaction with the outer parts of the disc leads to a larger Lindblad torque than that 
with the inner parts (Ward 1986, 1997). Therefore, the net Lindblad torque exerted by the planet causes it to lose angular 
momentum and to move inward relative to the gas (type I migration). Note that there is also a torque exerted by the planet 
on the material which corotates with the perturbation. Its sign depends on the gradient of vortensity at corotation (GT79). 
However, it is usually found to be less significant than the net Lindblad torque. The drift timescale for a 1-10 planet 
undergoing type I migration in a standard gaseous disc is typically between 10 4 and 10 6 years, shorter than the disc lifetime 
or estimated planetary formation timescales. 

When the mass of the planet becomes large enough, the perturbation becomes nonlinear. Feedback torques from the disc 
on the planet, which have to be taken into account in this regime, then opposes the motion of the planet and stops it altogether 
when the planet is massive enough (Ward 1997). However, density waves still transport angular momentum outward. If they 
are dissipated locally, the angular momentum they transport is deposited in the disc in the vicinity of the planet and a gap 
may be cleared out (Goldreich & Tremaine 1980; Lin & Papaloizou 1979, 1993 and references therein). The planet is then 
locked into the angular momentum transport process of the disc, and migrates inward at a rate controlled by the disc viscous 
timescale (type II migration). Here again, the drift timescale is rather short, being in the range 10 3 -10 5 years for standard 
parameters. 

Type II migration, which occurs for planets with masses at least comparable to that of Jupiter, can be avoided if after 
the planet forms there is not enough mass left in the disc to absorb its angular momentum. However, type I migration appears 
inevitable as there has to be enough gas in the disc to push a forming core inward if this core is to accrete a massive envelope 
to become a giant planet at some point. Different scenarii for stopping a planet undergoing inward migration have been 
considered (see, e.g., Terquem, Papaloizou & Nelson 2000 and references therein), but none of them can satisfactorily explain 
the presence of extrasolar planets with semi-major axes from a few 0.01 au all the way up to several au. 

Recently, Papaloizou (2002) has shown that type I migration reverses for reasonable disc models once the eccentricity 
of the orbit becomes comparable to the disc aspect ratio. This is because in that case the planet spends more time near 
apocentre, where it is being speeded up, than near pericentre, where it is being slowed down. Although the interaction with 
the disc tends to circularize the orbit of the planet (Goldreich & Tremaine 1980), significant eccentricity may be maintained 
by gravitational interactions between different planets forming simultaneously (Papaloizou & Larwood 2000). 

Here we investigate the effect of a magnetic field on planet migration in the linear regime. We consider a planet on a 
circular orbit and, to keep the problem tractable, restrict ourselves to the case of a purely toroidal field. This tends to be the 
dominant component in discs as it is produced by the shearing of radial field lines. As described above, the torque exerted by 
a planet on a disc depends mainly on the location of the radii where the perturbation is in resonance with the free oscillations 
of the disc. We therefore expect a magnetic field to modify the tidal torque, as it introduces more degrees of freedom in the 
disc. 

The plan of the paper is as follows. In sections 2, 3 and 4 we give the basic equations, describe the equilibrium disc 
model and give the expression of the perturbing potential acting on the disc. In section 5 we study the disc response to this 
perturbing potential. We first linearize the equations and derive the second order differential equation describing the disc 
response (§ 5.1). We then carry on a WKB analysis (§ 5.2), which shows that magnetosonic waves propagate outside the 
outermost turning points. In § 5.3 we study the disc response at the locations where it is singular. We call these radii magnetic 
resonances. There are two such resonances, located on each side of the planet's orbit and within the Lindblad resonances. In 
§ 5.4 we calculate all the turning points associated with the disc response. On each side of the planet there are two or three of 
them, depending on the value of the azimuthal number m. One of these turning points coincide with the Lindblad resonance. 
In contrast to the nonmagnetic case, waves are found to propagate in a restricted region inside the outermost turning points, 
around the magnetic resonances. We then proceed to give an expression for the tidal torque in general and in the vicinity 
of the magnetic resonances in particular in section 6. Numerical calculations of the torque exerted by the planet on the disc 
are presented in section 7. We show that the torque exerted in the vicinity of the magnetic resonances tends to dominate the 
disc response when the magnetic field is large enough. This torque, like the Lindblad torque, is negative inside the planet's 
orbit and positive outside the orbit. Therefore, if fi = c 2 /v\ increases fast enough with radius, the outer magnetic resonance 
becomes less important (it disappears altogether when there is no magnetic field outside the planet's orbit) and the total 
torque becomes negative, dominated by the inner magnetic resonance. This corresponds to a positive torque on the planet, 
which leads to outward migration. Finally, we discuss our results in section 8. 



2 BASIC EQUATIONS 

The disc is described by the equation of motion: 
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p(^ + (v- V)v) =-VP + F- /9 V*, (l) 

the equation of continuity: 

g! + V ■ (pv) = 0, (2) 
and the induction equation in the ideal MHD approximation: 

J = Vx(vxB), (3) 
where 

F=— (VxB)xB (4) 

Mo 

is the Lorentz force per unit volume, P the pressure, p the mass density, v the flow velocity, 9 the gravitational potential and 
B the magnetic field (fio is the permeability of vacuum). SI units are used throughout the paper. Here we neglect self gravity 
so that the gravitational potential * is assumed to be due to a central mass M* . 

We consider a thin disc, so that the above equations can be averaged over the disc thickness. To first order we neglect 
the 2-dependence of v and <]/, so that the vertically averaged equation of motion and continuity are, respectively: 



<9v 

W + (vV)v 



= -V(P) + (F) - EV#, (5) 
+ V • (Ev) = 0, (6) 



E 

dE 
dt 

where E is the surface mass density and the brackets denote averaging over the disc thickness. 
To close the system of equations, we adopt a barotropic equation of state: 

<P> = <P> (E) . (7) 

The sound speed c is then given by: 



3 EQUILIBRIUM DISC STRUCTURE 

We adopt a nonrotating cylindrical polar coordinate system (r,ip,z) with origin at the central mass. We denote (e r , e v ,e z ) 
the associated unit vectors. We suppose that at equilibrium the disc is axisymmetric and in rotation around a central mass, 
so that v = (0, rQ(r), 0), where O is the angular velocity. Furthermore, we assume that the equilibrium configuration contains 
only a toroidal magnetic field, i.e. B = (0, B(r, z), 0). The Lorentz force per unit volume is then: 

1 ( B °W ,-B%). (9) 



(in 1 r dr ' ' dz 

If we assume reflection with respect to the disc midplane and that B vanishes at the disc surface, this leads to: 

(10) 



4 PERTURBING POTENTIAL 

We consider a planet of mass M v <C M, on a circular orbit with radius r p and angular velocity Q p — \J GM, / r| . At the 
location (r, (p) in the disc, it exerts the gravitational potential: 

%{r, V ,t) = K(r,4-) - - — f M ; -ttj, (H) 

(Vq + r-p + r z — 2rr v cos <p) 

where <j> = ip — Q p t and we have introduced a softening length ro. We now expand *S>' p in a Fourier series with respect to the 
variable d>: 



<S>' p (rA) = yUr)cosm(f>, (12) 

m=0 

with 

2 r 

*UW = - / *pO.0) cosm4>d4>, 
n Jo 
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and 

*oW = - / %{r,4>)d4>. 
71 Jo 

Using the generalized Laplace coefficients defined as (Ward 1989; Korycansky & Pollack 93, hereafter KP93): 

,m , n 2 cos m<j> d(/> . . 

61/2(7) = - / YTT — -173 > (13) 

n Jo (q + P i - 27 cos 4>) 

where: 

■r ^ r 1 2 , . rl 

iir<r p : 7=—, p = 1, g = 1 + — , 

■r . r P * 2 * rfi 

if r>r p : 7=—, g = 1, p = 1 + — , 
we can rewrite ty' under the form: 



%n±o _ M p r p 



2 fl 2 M, r lf 



&r /2 (7) if>>r p , (14) 



=-s^» ifr< ^ (15) 

and the expression on the right-hand-side has to be divided by a factor two if m = 0. The subscript 'p' indicates that the 
quantities have to be evaluated at r = r p . 



5 THE DISC RESPONSE 

5.1 Linearization of the basic equations 

In this paper we consider small perturbations, so that the basic equations can be linearized. We denote Eulerian perturbations 
with a prime. Since the equilibrium is axisymmetric and steady, we can expand each of the perturbed quantities in Fourier 
series with respect to the variable <f> and solve separately for each value of m. The general problem may then be reduced 
to calculating the response of the disc to the real part of a complex potential of the form >I'J 7l (r)exp [i (rmp — cot)], where 
the amplitude ty' m is real and the frequency lj = mfl p , and summing up over m. We make all Eulerian fluid state variable 
perturbations complex by writing: 

oc 

X\r lV ,t) = ^X' m {r)S mv - ut \ (16) 

m=0 

where X is any state variable. The physical perturbations will be recovered by taking the real part of these complex quantities. 
We denote £ the Lagrangian displacement, and write its m-th Fourier component as £ m (r)exp [i (rrvp — tot)}. Since here we 
are only interested in the perturbations induced by the planet in the plane of the disc, we take £ z = 0. 

Linearization of the induction equation (3) , or the equivalent statement of the conservation of magnetic flux, leads to an 
expression for the perturbed magnetic field in the form (Chandrasekhar & Fermi 1953): 

B' = Vx(£xB). (17) 
The corresponding perturbation of the Lorentz force per unit volume is: 

F' = — |"(VxB') xB + (VxB) xB'l . (18) 
Mo 

Use of equations (17) and (18) with £ z = gives: 



IF' ) 



(B 2 ) 
Hor 



d £mr 
r 

dr 2 



, / 36l , \ df mr , ( 62 

+ (-- 1 j^T + U- m J — 



(19) 



{F ^ } - 1^7^^' (20) 

(FL) = 0, (21) 
where we have defined the dimensionless quantities: 



dlnr (B 2 ) dr \ dr 
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We now linearize the r- and (^-components of the equation of motion (5) and the equation of continuity (6). Using 
v' mr — ima^rnr, with a = O — O p , and equations (7) and (8), we get: 

w 



-m 2 a 2 £mr 



20v' 



<7 K~ „ 1 



-(F r ) + {FLr) 



<jW' 



Id v 
2 rL dr r 



(23) 

(24) 
(25) 



where we have found it convenient to use the radial component of the Lagrangian displacement £ mr and the azimuthal 
component of the Eulerian perturbed velocity v' mtp as variables. Here W' = E'c 2 /S is the linear perturbation of the enthalpy 
and k is the epicyclic frequency, which is defined by: 

2 20 d (r 2 0) 

K = j — ■ 

r ar 

We can further eliminate v' mip and W' m to get, after some tedious algebra, the following second-order differential equation 



for £„ 



d 2 i„ 



dr 2 



■ + 



Ai d£ 
r dr 



1 d<f' m 
c 2 dr 



with: 
A 2 = 

Ai = 



2c 2 



r 2 n 2 _ r 2 



K 20. 

20^ + ~ + 



ci + 1 + di + 



(26) 

(27) 
(28) 



-4o = 



20. 

^- I - d'i + d 2 + 1 - m 2 + 2cidi + — di 



~2iV + V- 1+Cl 



, 20 ^ 2 
1 + 

a r*o 



2_2 



bi 



20 
a 



di 
2 



+ 



r z er z — c 



k 20 
2?V + ~ 



1+Cl 



62 2 . 
+— - m + 



2 cr 2 / K 2 2Q , 

o + 1 



Cl 



fe 2 e 2 
4/3r 2 a 2 

2r 2 0cr 



2/3 



(29) 
(30) 

(31) 



r 2 cr 2 - c 2 y2fia 

where /? = c 2 /v\, with ua = \/ {B 2 ) / (/ioS) being the Alfven speed, and we have defined the dimensionless quantities: 

d _ dlnE d _ r^_(fZ c _ dine 
d In r ' £ dr 2 ' d In r 

The magnetic field comes in through the parameter j3 only. 

We see by multiplying all the coefficients of equation (26) by a 2 that there is no singularity in the coefficients of this 
equation where a = 0, i.e. at corotation, in contrast to the nonmagnetic case (GT79). When there is no magnetic field, the 
particles located at corotation can respond secularly to a forcing with frequency O = O p and therefore stay exactly in phase 
with the perturbation (they 'surf the tidal wave). In the presence of a toroidal field however, the tension of the field line 
to which the particles are attached provides a restoring force which introduces some inertia in the response of the particles, 
destroying the resonance. 

From equations (20), (24) and (25), we can write W' m and v' mlp as functions of £ mr . We give these expressions as they 
will be needed below: 



W' = 
with: 



dr r 



C = _^_[l n( Er C T)] + ^- 



b lC 2 



dlnr 



20a 2f3r 2 a 2 



(32) 



(33) 



and: 
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c d£ 



ra dr 



rcrc 



with: 



dlnr 



(34) 



(35) 



5.2 WKB dispersion relation 

We assume that the first-order derivative term in equation (26) can be neglected (this is verified a posteriori) and we define 
x = r/r v — 1, where r p is the corotation radius, i.e., for a circular orbit, the planet orbital radius. Then the homogeneous 
equation associated with equation (26) can be rewritten as: 

, . 1 ( , <? \\ d 2 ( £mr-\ A» £ mr , , 

1 + p{ 1 -^)\d^{—) + ^TW— = - (36) 

We assume \x\ <JC 1 and c <g; r \a\, which means \x\ >• H/r, where H ~ c/fi is the disc semi-thickness, and /3 ~ 1. With 1 + x 
being typically on the order of unity, the WKB approximation can then be used if \Ao\ 3> 1- We have: 

2 2 

AT °~ I 2 2 / 2\ 

~ —p - [m - k /a ) , 

providing 1 — k 2 / (m 2 cr 2 ) 3> c 2 / (r 2 a 2 ), i.e. a; 2 3> H 2 /r 2 + 1/m 2 . When this condition is satisfied then |.Ao| 3> 1 and we may 
look for homogeneous free-wave solutions of the form £ mr (r) = A(r)e lkr , where fcr ^> 1 and A is a slowly varying amplitude. 
Since _4i is of order (H/r) 2 x~ 3 , the first-order derivative term can be neglected compared to the second-order derivative 
term in equation (26) if x 3 S> (i?/r) 2 (fcr) _1 . When this is satisfied, the assumption made in this subsection is validated. Then 
equation (36) yields: 



-r-^i (37) 

which is the dispersion relation for magnetosonic waves (e.g., Tagger et al. 1990) propagating outside the Lindblad resonances 
(located at r = where m 2 a 2 — k 2 = 0) with the group velocity: 

k (c 2 + v\) 

v g = i (38) 

mo 

In the absence of magnetic field, these waves reduce to density waves (GT79). 

These solutions describe the response of the disc in the parts where the perturbing potential is negligible, i.e. outside the 
Lindblad resonances and far enough from them. 

We now consider the response of the disc inside the Lindblad resonances. 



5.3 Magnetic resonances 

The differential equation (26) is singular at the radii tm where A2 vanishes, i.e. where 1 + [l — c 2 / (?" 2 o" 2 )] / P = 0. We 
call these locations magnetic resonances. There are two such radii, the inner magnetic resonance r IM R < r p and the outer 
magnetic resonance tomr > r p . At these locations: 

m 2 (n-n p ) 2 = jfff A 2V (39) 
r i (v z A + c z ) 

i.e. the frequency of the perturbation in a frame rotating with the fluid matches that of a slow MHD wave propagating along 
the field line. 

Using Q, — Qk [l + O (// 2 /r 2 )] , where Ox is the Keplerian angular velocity, and H/r <Clwe get, to first order in H/r: 

1 1 2// . , „. 

\^-r P \ = ^—, (40) 

where H and f3 are evaluated in the vicinity of r p . We suppose here that these quantities vary on a scale large compared to 
H. As the field becomes weaker (/3 — » 00), the magnetic resonances converge toward the corotation radius. Note that since 
the effective Lindblad resonances are all located beyond 2H/3 from corotation (see section 5.4), they are beyond the magnetic 
resonances. 

We define the new variable x = (r — tm) /tm, where tm is either the inner or outer resonance, and we do a local analysis of 
equation (26) around x = 0. For \x\ <C 1, we have A2 — (dAij 'dr) r =r M (r — tm)- Here we assume that h = c/ (rfi) ~ H/r <g; 1, 
the parameters 61, 62, ci, di and e?2 are at most of order unity and (5 ~ 1. This amounts to assuming that only a varies in 
the expression (27) of Ai- For \x\ <JC 1, we then have: 
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A 2 ^ 2e 



(1 + /3m) 3/2 / dlnfi 



0. 



( dlnll \ x 
m \ d In r J rM h M ' 



where the subscript M indicates that the corresponding quantities are taken at r = rM, and e = — 1 or +1 depending on 
whether r M = r IM R or r MR- 

We now approximate Ai, Ao and So by their values at r — vm- This requires \x\ <JC hm- If this condition is not satisfied, 
Q — Q p cannot be replaced by its value at r = tm- To the first non zero order in Hm and allowing for m ~ h M x , we can then 
rewrite equation (26) under the form: 

where we have defined the dimensionless quantity £ mr = £ mr /rM and where A and 5 are given by: 



(44) 



3(l + /3 M ) 3/2 hAf V 1 + ^ M ^ M 

5 = 



3 (1 + /3m) /»m 



(1 + /3m) 1/2 \ r2Q2 9x J x=0 h M V r2f22 



e/3 M ( 1 d*^ \ 3 + (3 M ( *™ 



We next introduce the new variable u = 2\/Ax, so that equation (42) becomes: 

d £mr 1 d^mr z S . . 

+ r^-+^ = ^r- ( 45 ) 



9m 2 u du A 

Note that m is real outside the magnetic resonances, where Ax > 0, whereas it is imaginary inside the resonances, where 
Ax < 0. The solutions of this equation can be written as: 

Ur («) = J + ^C € Y (U) + C' e Mu), (46) 

where Jo and Yo are the Bessel functions of first and second kind, respectively, and C e and C' e are (complex) constants which 
a priori depend on /3m, Hm and m. The factor 7r/2 has been introduced to simplify the asymptotic expression of £ below. 

When u — » (i.e. |x| <g; /im), Jo(m) ~ 1 and Fo(m) ~ (2/ 7r ) lnu (Abramowitz & Stegun 1972). To deal with the singularity 
at u = 0, we use the Landau prescription, i.e. we displace the pole slightly from the real axis by replacing f2 — f2 p with 
Q. — £l p — iT, where V is a small positive real constant. This amounts to multiplying all the perturbed quantities by exp (Ft), 
and therefore to considering the response of the disc to a slowly increasing perturbation. This procedure can also be thought 
of as introducing a viscous force, or as shifting the resonances progressively with time. Note that such a shift occurs when the 
planet drifts through the disc, since the material in the resonances then changes with time. Shifting f2 — £l p below the real 
axis is equivalent to shifting x above this axis in equation (42) , since only — Q p was allowed to vary in the expression of A2 ■ 
Therefore, we replace x with x + ry, where 7 is a small positive real constant. For \x\ <C /im, the solutions of equation (45) 
can then be written as: 

imr ^ In [4A (x + i-y)] , (47) 
so that, near x = 0: 

Imr ~ C\ (\n \4A-y\ + i arctan ^ . (48) 

Note that there is a phase shift of tv on passing through x = 0. 

Since S/A has opposite sign at r — riMR and r = roMR (d^' m /dx changes sign on passing through corotation whereas 
does not) and is real, the real part of the solutions of equation (45) should have opposite sign at r = riMR and r = roMR 

whereas their imaginary parts should have the same sign. Therefore, from equation (48) we get C+i = — C*i, where the 

asterisk denotes the complex conjugate. 

Equations (32) and (34) can now be used to evaluate v' mlp and W' m for \x\ <C Km'- 



W 

y v m. 



a _ , ~~ . . , (49) 



r M n li 6 /3m (21 + 17)' 

"my ^ eC (1 +0M) 1/2 h M , 5Q s 

tm^m e 0m {x + i-y) 

where Q,m is f2 taken at r = rM- Finally, yields, for \x\ <JC Km'- 

_ mfe M ln[4A{x + il)l (51) 



tmVm (1 + 0m) 1/2 

We note that, in a nonmagnetized disc, the azimuthal perturbed velocity varies as In a; with x — \r — r p \ /r p — » near 
corotation (GT79), whereas here it varies as 1/x near the magnetic resonances. 
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5.4 Turning points 

At the locations where the solutions of equation (26) change from wave-like to non wave-like (turning points), the disc 
response varies on a large scale, and therefore it couples well to the perturbing potential. This results in a large contribution 
to the tidal torque (see GT79 and § 6). It is therefore important to calculate the location of these turning points. 

The WKB analysis above (see § 5.2) indicates that the differential equation (26) has turning points at the locations of 
the so-called nominal Lindblad resonances r — tl- Waves propagate outside these resonances, which are located on both 
sides of the corotation radius. In the nonmagnetic case, it is known that the turning points do no coincide with the nominal 
Lindblad resonances for values of m in excess of r/H (Goldreich & Tremaine 1980, Artymowicz 1993). If they did coincide, 
the turning points (also called effective Lindblad resonances) would converge toward r p as m increases. Instead, they converge 
toward a radius which is located a finite distance (equal to 2H/3) away from r p . Since the potential is more and more localized 
around r p at large m, coupling between the perturbation and the disc response is lost, and the contribution to the total torque 
becomes negligible. This is the torque cutoff first discussed by Goldreich & Tremaine (1980). 

The WKB analysis above indicates that for values of m up to about r/H, the location of the turning points does not 
depend on whether there is a magnetic field or not. In both cases they coincide with the nominal Lindblad resonances. 
However, for larger values of m, there are differences between a magnetized and a nonmagnetized disc, as the analysis below 
shows. We now calculate the location of the turning points. 



5.4-1 Elimination of the first-order derivative term in equation (26): 

There are no singularities in the coefficients of equation (26) at r = t\l. This is similar to the nonmagnetic case, where only the 
second-order differential equation for W' m has (removable) singularities there (Narayan, Goldreich & Goodman 1987). All the 
physical quantities are regular at these locations. To identify the turning points, we eliminate the first-order derivative term 
in this equation. We define x = (r — rh)/rL, tl being either the inner or outer nominal Lindblad resonance, and introduce 
the new dependent variable: 

y = £mr exp J dx^j , (52) 
so that the homogeneous equation associated with equation (26) can be written under the exact form: 

0+JCy = O, (53) 
where: 

*-£-i(£)'-5=(£)- 

We now expand K. around x = 0. At r = rL, k 2 = m 2 a 2 . We neglect here the departure from Keplerian rotation (which 
is second order in H/r) so that k = Q. Then a = e£l(rL)X/m, with X = 1 + 3emx/2 and e = +1 or —1 depending on whether 
we consider the outer or inner resonance, respectively. We assume (5 and h = c/(rSl) ~ H/r constant. Note that for values of 
m ~ 1/h and for x ~ h, we cannot approximate a by its value at tl- After some algebra, we find, to the first non zero order: 

— r- A 2 -l 6 27 m 2 h 2 

A\ ~ m 2 h 2 X 2 -m 2 h 2 4 (j2-m 2 h 2 ) 2 



HO 



m 



X 



X 2 - 1 „ 6 I f m 2 h 



m?h? X 2 -m 2 h 2 p\ X 



9m h / 3 , mh 

2X 2 (X 2 - m?h?) \2 2X 2 ' » 

where we have assumed m^> 1, since these are the values of interest for the tidal torque (see GT79). This expression is valid 
even for m larger than or very large compared to 1/h, and whatever the values of bi, 62, ci, d\ and <fe compatible with (3 and 
h being constant and provided they are not large compared to unity. Note that the expression above can also be seen as an 
expansion around the corotation radius, which corresponds to X = 0. 



5.4.2 Calculation of the turning points: 

The turning points are the radii where K, — 0. We solve this equation numerically and display the results in figure 1, where we 
plot (r — r p )/H (where H = c/fi) as a function of m for the turning points located beyond corotation, and for f3 = 0.1, 1 and 
10 and H/r = 0.03. The outer magnetic resonance is also represented. To each of these outer turning points and resonance 
corresponds an inner turning point or resonance located at the same distance from corotation to first order in H/r or 1/ra. 
Similar plots are obtained for other values of H/r. 

Artymowicz (1993) has calculated that, in a nonmagnetized disc, the location of the turning points is given by (his 
eq. [39]): 



Stopping migration by a toroidal magnetic field 9 



so that, for m — > oo, the distance to corotation is 2H/3. 

In a magnetized, we identify up to three turning points, that we label Rl, R2 and R3: 

(i) The outermost turning point Rl coincides with the nominal Lindblad resonance for small values of m. When m — ► oo, its 
distance to corotation is 2H/3 for j3 > 1 (like in the nonmagnetic case) and 2H / (3^/0) for (3 < 1, i.e. 2H' /3 with H' = va/Q- 

(ii) The intermediate turning point R2 is located at a distance from corotation equal to 2H/(3\f(3) for (3 > 1 and 2H/3 for 
(3< 1. 

(iii) The innermost turning point R3 exists only above some critical value of m, m cr i t , which depends on (3 and i//r. 
Interpolation of the numerical results gives m cr i t ~ 1.5v^ x r/H. When m — > oo, the distance of this turning point to 
corotation is 2H/(3\/P + 1), i.e. it coincides with the magnetic resonance. For either (3 <C 1 or /3 S> 1, R2 and R3 merge at 
large m. 

In figure 1, the shaded areas indicate where the waves are evanescent (i.e. K, < 0). In contrast to the nonmagnetic case, 
there is a region inside the outermost turning points where waves propagate. This region moves toward the corotation radius 
and decreases in size as (3 increases. It disappears altogether for infinite f3. Note that the magnetic resonances are contained 
within this region, so that the singular slow modes excited at the magnetic resonances can propagate. 

Since waves do not propagate around corotation for m > m cr i t , we expect a torque cutoff at large m like in the nonmagnetic 

case. 

We can check whether waves propagate or not around corotation by taking the limit X — > (i.e. r — > r p ) in the 
expression (55) above. In the nonmagnetic case (f3 — > oo): 

m 2 „ , 7 



£ — > -1 



A\ 4m 2 h 2 L ' 

Since this is negative, the waves are evanescent around corotation in that case. When f3 is finite: 

m 2 K t (m 2 hl) 2 ( 1 | 9 



X 4 V ^ <±™?h L 2 

This is positive when m < m cr it with m cr it = 3\ff3/(2hL), so that waves propagate around corotation for these values of m. 
Note that the value of m cr i t found by interpolation of the numerical results agrees perfectly with the analytical expression. 



6 TIDAL TORQUE 

6.1 Expression for the torque 

In the linear regime, the net torque exerted by the planet on the disc between the radii r\ and r-z is: 

T(n,ra) = ^ T m (n,r 2 ), 

with: 

T m (n,r 2 ) = ~J 2 j o Re S + Slr^"" 1 ) rxRc v(*L ™ (r!1,t) ) rd^dr. 



(56) 



Because of the ^-periodicity, the first-order term is zero and the torque has only a 2-component which can be written as: 

r 2 dr 

T m (n,r 2 ) = 2n ^ rdr (57) 



where, following KP93, we have defined the "torque density": 

dT m _ m /w* lT ,/ \ 
— = T Im(S m * m ). 

The asterisk denotes the complex conjugate. 

In a nonmagnetized disc, another useful quantity is the angular momentum flux: 



F m (r) = Er : 



f 

Jo 



Re 



im{ip — Q. p t^ 



Re 



; im(ip-n p t) 
u m,r 



which can also be written under the form: 
F m (r) = 7rSr 2 Re (y' mT v^) . 



(58) 



(59) 



(60) 
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In a nonmagnetized disc, inside the Lindblad resonances and away from corotation, the response of the disc is non 
wave-like and in phase with the perturbation, so that E^^^ is real. Outside the Lindblad resonances, the wavelength of the 
perturbation is small compared to the scale on which the perturbing potential varies so that rY/^'m integrated over radius 
is small. Therefore, most of the contribution to the torque comes from the resonances (Lynden-Bell & Kalnajs 1972, GT79). 
There are two contributions to the total torque. These are the Lindblad torque, which is equal to the angular momentum 
advected through the disc, and the corotation torque, which accumulates at the corotation resonance (GT79) and is associated 
with a discontinuity of the angular momentum flux there. 

In the presence of a magnetic field, it is not clear it is possible to calculate the torque separately at the different turning 
points, because the scale on which the disc response varies there may become larger than the distance between the turning 
points. In addition, the torque exerted around corotation for the values of m for which waves propagate there may be 
significant, as the scale on which the perturbing potential varies and that on which the disc response varies may be similar. 
The total torque will then be calculated numerically. 

However, an asymptotic form of the torque density at the magnetic resonances can be obtained, and this is done in the 
next subsection. 

6.2 Torque at the magnetic resonances 



The torque density in the vicinity of a magnetic resonance is, from equation (58) and with T,' m = TjW^/c 2 given by equa- 
tion (49): 



To get the torque exerted over a magnetic resonance, we integrate the torque density over the width Ar of the resonance, 
i.e. from tm — Ar/2 to tm + Ar/2. Note that Ax = Ar/rM ~ 7- Equation (49) has been used to derive the above expression 
for the torque density. Therefore, only if Aa;, and hence 7, is small compared to Km can this expression be used to calculate 
the torque over the resonance. When this condition is satisfied, the torque can be written as: 



with arctan [Ax/(27)] ~ 7r/2 if Aa: is a few times 7. 

Here we have assumed that the resonance were symmetrical around vm, so that the term proportional to Im(C*) in 
equation (62) gives no contribution. Since the integral over the term proportional to Re(Ce) converges (does not depend on 
the width of the resonance), the torque we obtain can be considered as a point-like contribution from the resonance, similar 
to the corotation torque in the nonmagnetic case. However, the resonance may not be symmetrical, in which case the term 
proportional to Im (C*) gives a global contribution, which does depend on the extent of the resonance. We will see below that 
the torque around the magnetic resonances is indeed dominated by this global contribution, and not by the point-like torque. 

Since Re(C+i) = — Re(C_i), the torque derived above has opposite sign at the inner and outer magnetic resonances. 

To know the amplitude and the sign of we need to calculate Im(C e ). This cannot be done from the local analysis 
presented in section 5.3, as only the solutions of the homogeneous equation determine the behaviour of the disc response 
at the magnetic resonances. The solutions obtained around the magnetic resonances are valid only for \r — tm\ <S H and 
therefore cannot be matched to the WKB solutions calculated in section 5.2. The constants C* e will then be calculated by 
matching the numerical and analytical solutions, and the asymptotic form of the torque density given by equation (61) will 
be used to obtain a precise value of the torque at the magnetic resonances. 

7 NUMERICAL CALCULATIONS 

In order to calculate the tidal torque exerted by the planet on the disc, we have solved the second-order linearized differential 
equation (26) numerically. 

7.1 Numerical scheme 

The numerical scheme we have used is similar to that presented by KP93. 

The coefficients of the Fourier decomposition of the potential (see eq. [12]) are calculated using the recurrence formulae 
between the generalized Laplace coefficients (see the appendix of KP93, which unfortunately contains several typos). 

Equation (26) is then integrated using the fifth-order Runge-Kutta procedure with adaptative stepsize control given by 
Press et al. (1992) from the disc inner radius r in to its outer radius r out . Note that KP93 started the integration at the 




(61) 



(62) 




(63) 
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corotation resonance and shot toward the boundaries. This procedure would not be appropriate here as in the presence of 
a magnetic field there is a resonance on each side of corotation. For a given value of m, we first compute two independent 
solutions £ mr ,i and £mr,2 of the homogeneous equation, and one particular solution £ mr , P of the inhomogeneous equation 
starting with random initial values. The general solution is then given by £ mr = £mr, P + ci£ mr ,i + C2^mr,2, where ci and C2 
are constants which depend on the boundary conditions. 

To calculate Ci and C2, we use the WKB solution of the equation, which should be valid at r in and r out if the boundaries are 
far enough from corotation. Then, at these locations, d^mr/dr — ik£ mr , where k is given by equation (37). Since the physical 
solutions are outgoing waves, we take k > 0. The group velocity is then positive (negative) outside (inside) corotation. These 
two boundary conditions yield two algebraic equations for ci and C2 (see eq. [24] of KP93) which can be solved to find these 
coefficients. 

To get better numerical accuracy, we repeat the procedure taking for the starting value of the particular solution £ mr , P 
the value of £ mr calculated at r — r in with the above coefficients ci and c 2 . This leads to £ mr , P converging toward the full 
solution £ mr . When recalculating ci and C2 after one iteration we indeed verify that these two constants are very close to zero. 
This procedure prevents large cancellation between £ mr , p and the general solutions of the homogeneous equation. 

Note that KP93 found the WKB approximation not to be accurate enough for the boundary conditions, and they 
accordingly took into account the amplitude as well as the phase variation of the WKB solution. We found that we could 
reproduce their results without having to go beyond the WKB approximation. 

Once £ mr is known, W' m and v' m(fi can be calculated from equations (32) and (34), v' mr — ima^ mr , and the torque density 
dTm/dA is given by equation (58). To calculate the total torque T m exerted by the planet on the disc, we integrate dT m /dA 
over the disc surface (see eq. [57]). 

7.2 Numerical results 

The planet potential softening length is set to r = 10~ 4 r p . Note that in the presence of a magnetic field the results are 
insensitive to the value of ro provided it is small enough, since there is no singularity of the homogeneous equation at r — r p , 
where the potential is singular. The Landau parameter is set to 7 = 1CP 6 . As expected, the total torque does not depend on 
the value of 7 provided it is small enough. 

For m < 5, we integrate the equation from r in — 0.2r p to r out = 5r p , whereas, for m > 6, we limit the range of integration 
from r in = 0.5r p to r out = 1.5r p . 

We assume that E, c and r 2 (B 2 ) vary like power laws of r. According to equations (22) and (31), we then have S oc r dl , 
c oc r ci and r 2 {B 2 } oc r bl . Then d 2 = - 1) and b 2 = (61 - 2) (61 - 1). Note that H/r ~ c/ (rfi) oc r c i+ B and 

= <? /v\ oc r 2c i+ d i _6 i+ 2 . We comment that the values of the parameters we chose below may require a departure from 
Keplerian rotation that is second order in H/r to ensure hydrostatic equilibrium. This is here neglected, as in Ward (1986, 
1997). 

To simplify the discussion, we define the dimensionless quantities W' m = W m j (r p Q p ^j, f m = T m / (T, p r p Q 2 ^, T m ,D = 
2n(T m /dA)/ (£ p r^) and F m = F m / (T, p r p Q. 2 ^, where the subscript p indicates that the quantity is taken at r = r p . The 
numerical results presented in this section correspond to M p /M* = 1. Since ^' m /(r p Q 2 ) oc M p /M* (see section 4), then 
W' m oc M p /M, and f m oc (M p /M,) 2 . 

7.2.1 CaseB = 

By setting B = 0, we have checked that we recover the KP93 results. In that case, the corotation torque is non zero if 
di =fc —1.5 and corresponds to a discontinuity of the angular momentum flux at r = r p . In figure 2, we plot W' m , the angular 
momentum flux F m (r) at radius r and the torque T m (r p ,r) exerted by the planet on the disc between the radii r p and r for 
B = 0, di = 0, c/ (r p Q p ) — 0.03 (ci = 0) and m = 10. For comparison, we also plot the angular momentum flux corresponding 
to di = —1.5 in the corotation region. The parameters used here are the same as those used by KP93 in their figures 2, 5 
and 6, so that a direct comparison can be made. For di = 0, the dimensionless corotation torque is about 102, whereas the 
dimensionless total torque T m (r in , r out ) is about 22, in very good agreement with the results displayed by KP93. 

7.2.2 Case B 

We first consider the case d\ = 0, c/ (r p fl p ) — 0.03 (ci = 0), 61 = 2, (3 = const = 1 and m = 10. In figures 3, 4 and 5 we 
plot W' m , Tm.,D and T m versus radius. Note that the rapid variation of Re(W^ l ) at r = r p is due to the rapid variation of the 
potential there, not to a singularity of the homogeneous differential equation. We have matched W m obtained numerically 
with the analytical expression (49) in order to compute the constants C e . We find that Re(C+i) > whereas Re(C_i) < 0. 
According to equation (63), the point-like torque exerted at the outer magnetic resonance is then negative whereas that 
exerted at the inner magnetic resonance is positive. This is in agreement with the curves shown in figure 4 around resonances. 
Figure 5 shows the total torque computed using the numerical form of the torque density. As a check, we have also calculated 
T m by using the analytical expression (61) for the torque density around the resonances (represented by the dashed lines in 
the lower panels of fig. 4). In general the two calculations are in good agreement with each other, although in some cases the 
discrepancy may reach about 50% (the sign of the torque does not change however). 
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In figure 6, we plot the total torque f m exerted by the planet on the disc versus m for different values of di , c\ and 61 and 
for both 5 = and B 7^ 0. In the magnetic case, the value of (5 at r p is taken to be 1. The cumulative torque T is obtained 
by summing up over m. We get an estimate by assuming f m is constant in the intervals 20 < m < 24, 25 < m < 29 etc... and 
f m = for m > 85. 

The values computed for B — are in very good agreement with KP93. In this case, the cumulative torque exerted by 
the planet on the disc is positive, which means that the planet loses angular momentum and migrates inward. When the 
magnetic field is non zero though, T may become negative, with a subsequent outward planet migration. This occurs when 
61 is negative enough, i.e. when /3 increases fast enough with radius or, equivalently, when the magnetic field decreases fast 
enough with radius. For instance, when di — ci = 0, T < for 61 = — 1, i.e. (3 oc r 3 or, equivalently, (B 2 ) oc r -3 , and for 
61 = 0, i.e. [3 oc r 2 or, equivalently, (B 2 ) oc r~ 2 . Note that this variation needs only be local. When &i goes from positive to 
negative values, the cumulative torque, initially positive, then becomes negative passing through zero. Therefore migration is 
slowed down and then reversed. 

In figure 7, we plot the one-sided torque exterted by the planet on the disc (i.e. the torque exterted inside and outside its 
orbit) versus m for di = and c/(r p Q p ) = 0.03, i.e. ci = 0. These plots show that the one-sided torque has a larger amplitude 
when a magnetic field is present. 

In figure 8 we plot the torque exerted between r in and r versus r/r p for m = 10, di = ci = 0, /3(r p ) = 1 and for 61 = 7, 2 
and —8. This shows that the total torque becomes negative when 61 is negative due to the contribution from the region inside 
the Lindblad resonances. 

In table 1 we have listed the value of the torque exerted on different radii intervals in the disc for m = 10, di = 0, 
c/(r p Q p ) = 0.03 (ci = 0), /3(r p ) = 1 and different values of 61. The torque exerted between the inner Lindblad resonance tilr 
and corotation r p is negative, whereas that exerted between corotation and the outer Lindblad resonance volr is positive. 
(Note that for m = 10 the Lindblad resonances coincide with the outermost turning points). We have found above that the 
point-like torque exerted at the inner magnetic resonance is positive whereas that exerted at the outer resonance is negative. 
Therefore, the torque is not dominated by the point-like contribution at the magnetic resonances, but is determined by the 
whole region around these resonances, where waves propagate for m < m cr ; t (see also the discussion in section 6.2). When bi 
decreases, j3 increases outward faster, so that the magnetic field becomes less and less important outside the planet's orbit. 
Therefore the torque exerted between r p and tolr, which is positive and tends to be negligible in the absence of a magnetic 
field, becomes less and less important. In contrast, the magnetic field becomes more important inside the orbit, where it gives 
rise to a negative torque. The net effect of decreasing 61 is then to decrease the total torque. If the magnetic field is strong 
enough, the total torque exerted on the disc is dominated by the region inside the Lindblad resonances (which is closer to the 
planet than the regions outside these resonances), and it can then become negative for small enough values of 61. 

Figure 9 shows the total torque exerted by the planet on the disc versus m for d\ = 0, c/(r p Q p ) = 0.03 (ci = 0), bi = — 1, 
which correspond to f3 oc r 3 , and different values of (3{r p ). As we can see, the cumulative torque T is significantly decreased 
for values of /3(r p ) as large as 10 2 and becomes negative for values of /3(r p ) between 10 2 and 10. 

The cumulative torque corresponding to f3(r p ) = 10 2 becomes negative if, for the same parameters as above, di is taken 
to be at least as large as unity, i.e. if E increases with radius at least as fast as linearly. This is because when di increases 
(5 increases faster with radius, which favours the negative inner torque exerted near corotation. This effect can be seen in 
figure 6 by comparing the upper and middle panels, which show that the total torque corresponding to di = —1.5 is larger 
than that corresponding to di = for 61 = — 1. Note that di > tends to favour the outer Lindblad resonance compared to 
the inner one, which results in a larger Lindblad torque (Ward 1986, 1997). The fact that the total torque actually decreases 
when di increases means that the variations of E affect more significantly the torque exerted around the magnetic resonances 
than the Lindblad torque. (Note that the variations of the torque displayed in figure 6 when d\ varies in the nonmagnetic 
case are mainly due to changes in the corotation torque). 

An increase of ci (corresponding to c or, equivalently, H/r, increasing faster with radius) also results in a smaller total 
torque, as can be seen by comparing the upper and lower panels of figure 6. When ci gets larger, /3 increases faster with 
radius, which again lowers the torque exerted around the magnetic resonances. Also the Lindblad torque is reduced when ci 
gets larger (Ward 1986, 1997). 

Figure 10 shows the effect of varying the value of H/r = c/(rQ) at r = r p for m = 10, d\ = 0, ci = 0, f3(r p ) = 1 and 
bi = — 1 and 3. The total torque exerted on the disc tends to zero as H/r increases. This is because the resonances are further 
away from the planet when H/r is larger, so that the coupling between the perturbation and the tidal potential is weaker. 
For 61 = — 1, the torque is always negative and decreases as H/r decreases. For b\ = 3, the torque is positive and maximum 
for some value of H/r (~ 0.03). It stays positive (and becomes smaller) as H/r increases, whereas it becomes negative for 
smaller values of H/r. This is because even though j3 increases with radius in that case, the region between the inner Lindblad 
resonance and corotation contributes to a larger magnitude of the torque than the region between corotation and the outer 
Lindblad resonances. This is also the case for H/r = 0.03 and 61 = 1, as shown in table 1. 

Note that what determines the sign of the torque is the gradient of f3, not of the magnetic field itself, as only (3 enters 
into equation (26). 
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8 DISCUSSION AND CONCLUSION 

We have calculated the linear torque exerted by a planet on a circular orbit on a disc containing a toroidal magnetic field. 
In contrast to the nonmagnetic case, there is no singularity at the corotation radius, where the frequency of the perturbation 
matches the orbital frequency. However, there are two new singularities on both sides of corotation, where the frequency of 
the perturbation in a frame rotating with the fluid matches that of a slow MHD wave propagating along the field line. These 
so-called magnetic resonances are closer to the planet than the Lindblad resonances. In addition, on each side of the planet, 
there are two turning points located beyond the magnetic resonance. One of them coincides with the Lindblad resonance. For 
values of m larger than some critical value m cr it, there is a third turning point between the magnetic resonance and corotation, 
on each side of the planet. Like in the nonmagnetic case, waves propagate outside the outermost turning points. But here they 
also propagate inside the intermediate turning points when m < m cr i t , or between the intermediate and innermost turning 
points when m > m cr i t (see figure 1). The singular modes excited at the magnetic resonances can therefore propagate. 

There is a significant torque exerted on the region of the disc inside the outermost turning points (which coincide with the 
Lindblad resonances for the values of m of interest). Like the Lindblad torque, this torque is negative inside the planet's orbit 
and positive outside its orbit. The whole region around corotation, not just a narrow zone around the magnetic resonances, 
contribute to this torque. In other words, the magnetic resonances contribute to a global, not a point-like, torque. Since these 
resonances are closer to the planet than the Lindblad resonances, they couple more strongly to the tidal potential. Therefore, 
the torque exerted around the magnetic resonances dominate over the Lindblad torque if the magnetic field is large enough. If 
in addition (3 = c 2 jv\ increases fast enough with radius, the outer magnetic resonance becomes less important (it disappears 
altogether when there is no magnetic field outside the planet's orbit) and the total torque is then negative, dominated by the 
inner magnetic resonance. This corresponds to a positive torque on the planet, which leads to outward migration. 

The amount by which /3 has to increase outward for the total torque exerted on the disc to be negative depends mainly 
on the magnitude of (3. We have found that for (3 = 1 at corotation, the cumulative torque (obtained by summing up the 
contributions from all the values of m) exerted on the disc is negative when (3 increases at least as fast as r 2 . If (3 oc r 3 , the 
cumulative torque becomes negative for values of (3(r p ) between 10 2 and 10, whereas it is negative for f3(r p ) — 10 2 if (3 oc r 4 . 

The migration timescales that correspond to the torques calculated above are rather short. The orbital decay timescale 
of a planet of mass M p at radius r p is r = M p r p Q. p /\T\, where T is the cumulative torque exerted by the planet on the disc. 
Remembering that T oc (M p /M*) 2 , we get: 

r(yr) = 4.3 x 10 9 ( *±) " ( —^—) " ( J>-) ^ f 
yy ' V M ®/ V 100 S cm / VlauJ 

In a standard disc model, E ~ 100-10 3 g cm~ 2 at 1 au (see, for instance, Papaloizou & Terquem 1999). Therefore, r ~ 10 5 - 
10 6 yr for a one earth mass planet at 1 au in a nonmagnetic disc, as T ~ 10 3 in that case (see fig. 6). This is in agreement 
with Ward (1986, 1997). In a magnetic disc, \T\ may become larger, leading to an even shorter migration timescale (note 
that given the limited accuracy of the numerical scheme used, as indicated in section 7.2.2, only orders of magnitude for the 
migration timescale are obtained here). However, it is important to keep in mind that these timescales are local. Once the 
planet migrates outward out of the region where (3 increases with radius, it may enter a region where f3 behaves differently 
and then resume inward migration for instance. Such a situation would be expected to occur in a turbulent magnetized disc 
in which the large scale field structure changes sufficiently slowly. 

Unless the magnetic field is above equipartition, it is unstable by the magnetorotational instability (Balbus & Hawley 1991, 
1998 and references therein), the saturated nonlinear outcome of which is MHD turbulence. In a turbulent magnetized disc, 
the main component of the field is toroidal, as shown by numerical simulations (Hawley, Gammie & Balbus 1995; Brandenburg 
et al. 1995). Global simulations have also shown that the turbulence saturates at a level corresponding to (3 ~ 100 if there is 
no mean flux or to lower values of (3 if there is a mean flux (Hawley 2001; Steinacker & Papaloizou 2002). Also, substantial 
spatial inhomogeneities are created by radial variations of the Maxwell stress (Hawley 2001; Steinacker & Papaloizou 2002), 
so that the field may display significant gradients. On the basis of the work presented here, we are led to speculate that in 
such a disc the planet would suffer alternately inward and outward migration, or even no migration at all. It would then 
oscillates back and forth in some region of the disc, or undergo some kind of diffusive migration, either outward or inward 
depending on the field gradients encountered. Note however that in a turbulent disc the magnetic field may vary locally on 
timescales shorter than the timescales needed to establish the type of tidal response assumed in this paper. 

It has been pointed out that protoplanetary discs may be ionized enough for the magnetic field to couple to the matter 
only in their innermost and outermost parts (Gammie 1996; Fromang, Terquem & Balbus 2002). A planet forming at around 
one astronomical unit, where the ionization fraction is very low, would then migrate inward on the timescale calculated by 
Ward (1986, 1997) until it reaches smaller radii where the field is coupled to the matter. At this point, there would be a 
magnetic field inside the planet's orbit but not outside its orbit. The torque on the planet would then reverse, and outward 
migration would occur. However, as soon as the planet would re-enter the nonmagnetic region, inward migration would 
resume. Hence, the planet would stall at the border between the magnetic and nonmagnetic regions. 

Inward migration in a magnetized disc may then either be very significantly slowed down, occur only on limited scales, 
or not occur at all. The planet would then be able to grow to become a terrestrial planet or the core of a giant planet. 

When the planet becomes massive enough (about 10 earth masses), the interaction with the disc becomes nonlinear, with 
a gap being opened up around the planet's orbit. Because a rather strong torque is exerted in the vicinity of the planet in 
the presence of a magnetic field, a gap should open up more easily and be 'cleaner' in a turbulent magnetic disc than in a 



(64) 
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similar laminar viscous disc. This seems to be in agreement with the numerical simulations of a giant planet in a turbulent 
disc performed by Nelson & Papaloizou (2002). Once the gap is completely cleared out, the effect of the magnetic resonances 
disappear as they are inside the edges of the gap (which coincide with the Lindblad resonances corresponding tom~r/F). 
The magnetic field may however still have an effect on the disc-planet interaction because of the turbulence it produces 
(Nelson & Papaloizou 2002). 

We have focused here on a planet with a circular orbit. A planet on an eccentric orbit in a nonmagnetic disc usually 
suffers eccentricity damping, as in that case the corotation resonances, which damp the eccentricity, dominate over the Lindblad 
resonances, which excite it (Goldreich & Tremaine f980; Papaloizou, Nelson & Masset 200f). This picture may change once 
a magnetic field is introduced. The effect of a magnetic field on the eccentricity of a planet will be studied in another paper. 
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Table 1. Torque T m (n, T2) exerted by the planet on the disc between the radii r\ and T2 in units S p rpf!p for different values of b± and 
for m = 10, d\ = 0, c/(r p Q p ) = 0.03 (ci = 0) and /3(r p ) = 1. The radii r\ and T2 are either the disc inner radius ri n , the inner Lindblad 
resonance rj^n, the corotation radius r p , the outer Lindblad resonance tqlr or the disc outer radius r ou t- For comparison, the last line 
shows the torques in the nonmagnetic case. As b\ decreases, the torque exerted inside the Lindblad resonances becomes negative and 
dominant. 
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Figure 1. Distance of the outer turning points (thin solid lines) and magnetic resonance (thick bold lines) to corotation in units of 
H = c/Q versus m for /3 = 0.1 (upper plot), 1 (middle plot) and 10 (lower plot) and for H/r = 0.03. The turning points are labelled 
Rl, R2 and R3. The shaded areas indicate the regions where the waves are evanescent. The location of the effective outer Lindblad 
resonance in a nonmagnetizcd disc is also represented (dashed lines). The inner turning points and resonance arc obtained by reflection 
with respect to the a;-axis. 
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Figure 2. Results for B = 0, d\ = d In S/dln r = 0, c/ (r p f2 p ) = 0.03 (ci = 0) and m = 10. Upper left panel: Real part (solid line) and 
imaginary part (dashed line) of W^, versus r/r p . Lower left panel: Tidal torque T m (r p ,r) exerted by the planet on the disc between the 
radii r v and r versus \r — r p \/r p . The two curves correspond to r > r p (solid line) and r < r p (dashed line). Upper right panel: Angular 
momentum flux F m versus \r — r p \/r p . The two curves correspond to r > r p (solid line) and r < r p (dashed line). Lower right panel: F m 
versus r/r p around corotation for d\ = (solid line) and -1.5 (dotted line). 
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Figure 3. Results for B ^£ 0: W' m versus r/r p for di = 0, c/ (r p U p ) = 0.03 (ci = 0), 6i = 2, /3 = const = 1 and m = 10. Both the real 
part (solid line) and imaginary part (dashed line) are plotted. The two lower panels show a zoom on the magnetic resonances (left panel, 
inner resonance; right panel, outer resonance). 
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Figure 4. Torque density T m ^ry versus r/fp for the same parameters as in fig. 3. The two lower panels show a zoom on the magnetic 
resonances {left panel, inner resonance; right panel, outer resonance). The dashed lines represent the fit given by equation (61), where the 
constants C e are calculated by matching W m obtained numerically with the analytical expression (49). The point— like torque exerted at 
the inner magnetic resonance is positive whereas that exerted at the outer magnetic resonance is negative. 
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Figure 5. Torque f m exerted by the planet on the disc for the same parameters as in fig. 3. The upper panel shows the torque exerted 
between the radii Ti n and r versus r/r p . The lower panel shows the torque exerted between the radii r p and r versus \r — r p \/r p , for 
both r > r p (solid line) and r <r p (dashed line). Here the total torque on the disc is positive, like in the nonmagnetic case. 
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Figure 6. Total torque T m (ri n ,r ou t) versus m for d\ = {upper and lower panels), and -1.5 {middle panel) and for c/{r p Q p ) = 0.03, 
i.e. ci = {upper and middle panels) and c/(rf2) = 0.03, i.e. ci = —0.5 {lower panel). The open squares correspond to /3 = oo, i.e. 
B = 0. The filled symbols correspond to (3(r p ) = 1 and bi = —1 {triangles), {pentagons ), 2 {squares) and 4 {hexagons), as indicated 
by the label on the curves. In the upper panel, an estimate of the cumulative torque gives T = —3975, -1300, 4160 and 9148 for B/0 
and 6i = -1, 0, 2 and 4, respectively, and f = 1362 for B = 0. In the middle panel, f = -1768 for 6i = -1 and T = 2051 for B = 0. 
In the lower panel, T = —2680 for 6i = — 1 and T = 2383 for B = 0. A negative (positive) cumulative torque corresponds to the planet 
migrating outward (inward). 
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Figure 7. Upper panel: minus the torque exterted by the planet on the disc inside its orbit, -f m (ri„, r p ), versus m. Lower panel: torque 
extorted by the planet on the disc outside its orbit, T m (r p ,r ut), versus m. The parameters are the same as in the upper panel of figure 6: 
di = and c/(r p Q p ) = 0.03, i.e. c\ = 0. Like in figure 6, the open squares correspond to j3 = oo, i.e. B = 0, and the filled symbols 
correspond to (3(r p ) = 1 and bi = —1 (triangles), (pentagons ), 2 (squares) and 4 (hexagons). An estimate of the cumulative torque for 
to < 80 gives {f(r in ,r p );f(r p ,r out )] = (-39531; 35560), (-37727; 36521), (-35442; 39535), (-32856; 42161) for B ^ and bi = -1, 0, 2 
and 4, respectively, and (—28278; 29486) for B = 0. The one-sided torque has a larger amplitude when a magnetic field is present. 
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Figure 8. Torque f m (ri„,r) exerted between r^ n and r versus r/r p for m = 10, di = 0, c/(r p n p ) = 0.03 (ci = 0), /3(r p ) = 1 and b\=7 
(solid line), 2 (short-dashed line) and -8 (long-dashed line). As b\ decreases, the torque exerted inside the Lindblad resonances becomes 
increasingly dominant and negative. 
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Figure 9. Total torque f m (ri n ,r ut) versus m for di = and c/(r p Q p ) = 0.03 (ci = 0). The open squares correspond to /3 = oo, 
i.e. S = 0. The filled symbols correspond to bi = —1 and /3(r p ) = 0.1 (triangles), 1 (squares), 10 (pentagons) 10 2 (hexagons) and 10 3 
(septagons). An estimate of the cumulative torque gives T = 1362, 1129, 855, -508, -3975 and -5483 for 0(r p ) = oo, 10 3 , 10 2 , 10, 1 and 
0.1, respectively. This shows that T is significantly reduced for values of /3 as large as 10 2 . 
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Figure 10. Total torque T m (ri n ,r ou t) versus H p /r p = c p /(r p Q p ) (where the subscript 'p' indicates that the quantities are taken at 
r = r p ) for m = 10, d\ = 0, c\ = 0, I3(r p ) = 1 and b\ = —1 (solid line) and 3 (dashed line). The torque tends to zero as H/r increases. 



